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Abstract 

We predict spontaneous nematic order in an ensemble of active force generators with elastic inter- 
actions as a minimal model for early nematic alignment of short stress fibers in non-motile, adhered 
cells. Mean-field theory is formally equivalent to Maier-Saupe theory for a nematic liquid. However, 
the elastic interactions are long-ranged (and thus depend on cell shape and matrix elasticity) and 
originate in cell activity. Depending on the density of force generators, we find two regimes of cellular 
rigidity sensing for which orientational, nematic order of stress fibers depends on matrix rigidity either 
in a step-like manner or with a maximum at an optimal rigidity. 
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I. INTRODUCTION 



The actin cytoskeleton of living cells comprises semiflexible actin filaments and force- 
generating myosin molecular motors (as well as many regulatory proteins) [1]. In different 
cell types, the actin cytoskeleton can display quite distinct degrees of order: In many motile 
cell types, the bulk cytoskeleton is simply a cross-linked filament meshwork without any higher 
degree of ordering, while in striated muscle cells, actin and myosin assemble into crystal-like 
myofibril aggregates [2]. We are interested in an intermediate state of cytoskeletal order found 
in non-motile, adhered cells, where crosslinked bundles of parallelly aligned actin filaments 
form, so called stress fibers [3]. These stress fibers often align with the long axis of the cell, 
which establishes nematic order within the ensemble of stress fibers and results in polarized cell 
forces. Interestingly, experiments demonstrate a strong influence of substrate rigidity on the 
nematic ordering of stress fibers. In recent experiments with stem cells, the degree of nematic 
order of nascent stress fibers showed a non-monotonic dependence on substrate rigidity with 
a maximum at an optimal rigidity [4]. More generally, substrate rigidity was shown to be a 
determining factor for morphology and could even trigger cell fate decisions [5]. 

In this paper, we address the question of why nascent stress fibers align preferentially with the 
long axis of the cell in a substrate dependent manner. Our work provides a unified theoretical 
foundation for a previous study, which addressed this question in the framework of a linear 
response theory that treated the entire cell as a single elastic inclusion [4]. We predict more 
general phase behavior for nematic ordering as a function of substrate rigidity, cell shape 
anisotropy and contractile force strength. We focus on the early stages of symmetry breaking 
by considering a minimal model of force generators that is motivated by short, nascent stress 
fibers, which are coupled to the remaining, disordered cytoskeleton and which exert contractile 
forces. The resulting elastic deformations of the cytoskeleton provide a mean of long-range 
communication between distant parts of the cell [6] . We propose that such elastic interactions 
between nascent stress fibers drive their nematic alignment with respect to the long axis of the 
cell in a way that depends on cell shape and substrate rigidity. 

More precisely, we study the emergence of nematic order within an ensemble of active force 
generators. We employ a generic description of the force generators in terms of extended force 
dipoles, which are subject to mutual elastic interactions [7, 8]. The force generators are ho- 
mogenously distributed within the cellular domain for which we assume suitable boundary 
conditions. For simplicity, we model the disordered cytoskeleton as an isotropic, linear elas- 
tic material. While the assumption of isotropy might be well justified for the early stages of 
cytoskeletal symmetry breaking, the assumption of a linear elastic material is a strong sim- 
plification: The cytoskeleton is known to be viscoelastic on long time-scales and may exhibit 
non-linear elastic behavior for large stresses. Despite these limitations, our minimal model pro- 
vides qualitative insight into cytoskeletal polarization and reveals different regimes of rigidity 
sensing. 

We find that the cellular shape determines a macroscopic strain field due to the long-ranged 
nature of elastic interactions. It thus provides a cue for nematic ordering similar to the shape- 
dependent depolarization factors in electrostatics [9] . However, in contrast to the electrostatic 
case where the shape modifies the response to an external electric field, the elastic case of ac- 
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tive force dipoles shows shape-dependent nematic order even in the absence of external stresses 
[10]. Additionally, generic, isotropic hard-core repulsion between force generators can favor 
cooperative nematic alignment within the ensemble. Using a mean field approach inspired by 
the classical theory of orientational polarization of interacting dipoles [9], we derive the de- 
pendence of nematic ordering on macroscopic boundary conditions (such as domain shape and 
rigidity of a surrounding matrix) from pairwise elastic interactions. This dependence provides 
our prototypical cell with a mechanism to sense and respond to substrate stiffness [11]. We use 
a mean field approach to predict nematic order of nascent stress fibers as a function of physical 
parameters of the surrounding matrix. In addition to a regime of non-monotonic dependence 
of nematic order on matrix stiffness with a maximum at some optimal rigidity, which had been 
already found in [4], we find a second regime, which is characterized by a monotonic, step-like 
dependence. Experiments on cells with different aspect ratios observed qualitatively different 
polarization responses that match the two regimes of our theoretical model. Our approach thus 
generalizes and unifies previous work that either described active cellular responses by a phe- 
nomenological linear response theory [4], or considered elastic interactions between individual 
point force dipoles, but without reference to cell shape or matrix rigidity [12]. 

II. EXTENDED FORCE DIPOLES 

We consider a simple model of a mesoscopic force generator embedded in an elastic material, 
which we propose as a generic description of the tension forces exerted by a short stress fiber 
on the surrounding, disordered cytoskeleton. A general force generator located at x = exerts 
a field /j(x) of active forces on the cytoskeleton and induces a strain field i%(x). The force 
field /j(x) is characterized by its multipole moments Pi lt ... t i nt j = f cPxx^ ■ ■ •Xj n / J -(x) [8]. Note 
that the monopole moment vanishes, Pj = 0, since the force generator is not acted upon by 
external forces. Thus the dipole moment P^ dominates the far field strain. The far-field 
strain at distances |x| much larger than the size a of the force generator equals to leading 
order the strain field induced by a point force dipole with localized force dipole density 
Pij(x) = Pij 5(x) [8]. If the cell were infinite, 



where E c is the Young's modulus of the cell and Q is the angular part of the force dipole Green's 
function that depends on its Poisson ratio v c [8, 13]. Note that the strain field is long-ranged 
and decays as |x| -3 just like the electric field of an electric dipole, but the symmetry of this 
field expressed by Q is akin to that of an electric quadrupole. The field does not accurately 
reflect the strain field at small distances comparable to the size a of the force generator. In 
particular, u^ r diverges at x = 0. Here, we abstract from fine details of the strain field 
at small distances and propose to employ a as a cut-off distance by introducing a force dipole 
density Pij(x) that takes the constant value Pij/\B\ inside a spherical region B = {|x| < a} 
of volume \B\ = 4ira 3 /3, but is zero outside B. Note that the force dipole density Pij{x) is 
equivalent to a set of surface forces f\ = PijUj acting on the boundary of B, where rij is the 
surface normal of the region B. This is analogous to the situation in electrostatics, where a 
uniform electric dipole density induces the same electrical field as a set of surface charges [9]. 




(1) 
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Using Eshelby's classical theory of elastic inclusions, we find that the strain field uf- inside the 
spherical region B is uniform [14, 15]. Its purely dilational strain u 1 ^' = u™ k 5ij/3 is linear in the 
isotropic part p^ = ptk of the dipole density, while the pure shear strain u^' s = — 
(which does not involve any volume change) is linear in the anisotropic part p\j = p^ — 

u7j'° = uoPij/Ec and w™' s = a s p\jlE c . (2) 

The coefficients a = (1/3) (1 + i/ c )/(l - u c )(l - 2v c ) and a s = (2/15) (4 - 5i/ c )/(l - i/ c )(l + i/ c ) 
are positive for all physical values of the Poisson ratio — 1 < v c < 1/2 [14]. For incompressible 
matrices characterized by v c = 1/2, the coefficient vanishes. The strain field u°J l outside 
the region B is equivalent to a strain field induced by the fictitious point force dipole Pjj5(x) 
plus some fictitious point force octupole P^^x) located at the center of the ball [15]. We will 
not require the explicit form of the outer strain field; for later calculations, it is sufficient to 
note that its spherical mean f| x | =c d 2 x-u° ut ' fc over a shell concentric with B vanishes. The inner 
and outer strain field are related by a jump condition that involves the fictitious surfaces force 
/ mentioned above, u™ — w° ut (an) = 3aofirij/E c for |n| = 1 [15]. Previous work considered 
point force dipoles, which correspond to the limit a — » 0. Our use of extended force dipoles 
avoids singularities, in particular the elastic deformation energy associated with a force dipole 
is finite. Additionally, hard-core repulsion between different force dipoles is introduced in a 
natural way. The extended force dipole provides a prototypical example of a force generator 
and will be used throughout this paper. 

We now consider an ensemble of homogenously distributed, extended force dipoles of density 
p that reside inside an elastic domain Q. The dipoles have random positions x^ and dipole 
moments P^ . The individual dipole moments P^ = P^n^n^ are bipolar with director n^. 
Such a dipole moment equals the dipole moment of a pair of opposing point forces ±P / an fc 
separated by a vector 2an.k. The elastic domain Q and the force dipoles mimic the contractile 
activity of short, nascent stress fibers embedded into disordered cytoskeleton within an non- 
motile, adhered cell. We assume that the domain Q is surrounded by an infinitely extended, 
elastic matrix that represents the physical substrate to which the cell adheres, see fig. 1. The 
Young's modulus E m and the Poisson ratio v m of this matrix may differ from the respective 
values E c and v c of the cellular domain Q. To capture the effect of an anisotropic domain 
geometry, yet keep the calculations simple, we consider spheroidal domains with semi-axis 
a x = a y < a z and aspect ratio r = a z /a x . Thus the z-axis denotes the axis of revolution of the 
cellular domain. We treat the dipole strength P as a control parameter; experiments suggest 
that Po m ay increase with matrix stiffness [16]. 

We assume that the directors n k of the individual force dipoles are free to rotate (while 
their positions x^ are fixed for simplicity). In our particular case of cellular domain that has 
the shape of a prolate spheroid, we can show that polarization only occurs along the axis of 
revolution of the cellular domain. We can therefore characterize nematic order by a single, 
scalar order parameter 

S = (3/2)(cos 2 ^), - (1/2), cos0 fe = n M . (3) 
This order parameter S takes the maximal value 1 for perfect alignment with the z-axis and 
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FIG. 1: We present a minimal model for the elastic interactions between short, nascent stress fibers 
within a non-motile, adhered cell. The short stress fibers are modeled as an ensemble of force generators 
that reside inside an elastic cellular domain £1 with Young's modulus E c and Poisson ratio v c . The 
cellular domain in turn is embedded in an elastic matrix (with Young's modulus E m and Poisson ratio 
u m ) that mimics the substrate to which the cell adheres. 

vanishes for an isotropic director distribution. For later use, we also introduce the alignment 
parameters Sk = (3/2) cos 2 9k — (1/2) of all the individual dipoles. 

III. ELASTIC INTERACTIONS 

The elastic deformation energy due to the ensemble of extended force dipoles is H — 
J <i 3 x (E/M^)(Sfccr^' ) )/2 where aff is the elastic stress field induced by the fc-th dipole [13]. 
Using the fact that the force due to an ensemble of dipoles is proportional to the gradient of 
the dipole density (Vjcr^ = V»pj^), we rewrite H as [8] 



Here, Ukk represents a "self-energy" of the k-ih dipole [24], whereas Uki represents a pairwise 
elastic interaction energy which can be written in the form Uki ~ P^ QimjnPmn [7, 8]. It 
has been proposed that active force generators that are fueled by an external energy reservoir 
minimize the work needed to deform the matrix while maintaining a constant force magnitude 
(here Pq/o) [17] [25]. This minimization principle successfully describes the behavior of whole 
cells adhered to an elastic substrate in a phenomenological way. Assuming that cytoskeletal 
reorganization is governed by local mechanosensitive feedback mechanisms, we argue that this 
minimization principle also applies to sub-cellular cytoskeletal structures such as nascent stress 
fibers. This is a crucial assumption of our work and can be tested by comparing our theory to 
future experiments. 
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IV. MACROSCOPIC STRAIN FIELDS DEPEND ON BOUNDARY CONDITIONS 



We introduce the macroscopic strain field Uij that averages the microscopic strain field 

(k) 

Sfc-u^- on a length-scale that is much larger than both the dipole ball size a and the typical 
dipole-dipole distance p -1 / 3 . Due to the linear relationship between the strain field and the 
stress sources, this macroscopic strain field is equivalent to the strain field induced by the 
macroscopic force dipole density p^ = pPij where P^ = (P-p} is the ensemble average of 
the force dipole moments [26]. The force dipole ensemble generates an isotropic macroscopic 
dipole density p~°- = p^ij/^ = PP)<%/3 irrespective of nematic ordering. This "hydrostatic 
pressure" has no analogue in the electrostatic case of electric dipoles, which have polar as 
opposed to tensor, nematic symmetry. A net polarization of the dipole ensemble (here assumed 
in the z-direction) gives an additional anisotropic part pf ■ = p^ — = pPoS(5i Z 5j Z — <%/3) 
which is proportional to the nematic order parameter S. For domains of spheroid shape, the 
corresponding macroscopic strain is homogenous throughout the domain Q [14], similar to the 
situation in electrostatics [9]. This homogenous strain can be split into a dilation part 
v%j = Ukk^ij and a pure shear part u\j = — TZ?-. Only the shear exerts a torque on a force 
dipole and is therefore the quantity of interest to us. From symmetry considerations, we infer 
that the shear strain can be written in the form 



The coefficients h and g can again be computed using Eshelby's theory of elastic inclusions [14] ; 
they are depicted in fig. 1(a) as a function of matrix rigidity E m for different values of the domain 
aspect ratio r. The coefficient g/E c couples the ^-anisotropic part of the macroscopic force 
dipole density to the macroscopic shear strain and thus provides a macroscopic analogue of 
the shear compliance l/(2/x c ) = (1 + u c )/E c . The coefficient h is novel in the elastic case and 
does not exist in the simple electrostatic case; it couples a "hydrostatic pressure" p° to a shear 
strain; this shape field factor must vanish by symmetry in the case of a spherical domain, but 
may be non-zero for non-spherical domains due an anisotropic distribution of elastic restoring 
forces from the matrix. Below, we infer the qualitative behavior of h and g by discussing four 
important limit cases: (i) The case of a matrix that is much stiffer than the cellular domain 
with E m ^> E c corresponds to clamped boundary conditions. Hence, the macroscopic elastic 
strain vanishes within the cellular domain and g = h = 0. In this case, restoring forces from 
the matrix completely counterbalance the active cell forces, (ii) If the matrix is much softer 
than the cellular domain, E m <C E c , the elastic domain Q has essentially stress-free boundary 
conditions, and there are no restoring forces from the matrix. Hence, the macroscopic elastic 
stress Uij within Q equals p^. From pf^ = (jf ■ = 2p c TZfj, we conclude g = 2p c E c — 1 + v c and 
h = 0, i.e. there is no shape-dependence of the macroscopic strain field, (iii) If SI is a sphere 
and has the same elastic properties as the matrix, then w^- = a s pP^/E c is computed analogous 
to the inner strain field of an extended force dipole. Hence, g — a s , h = 0. (iv) Finally, in the 
limiting case of a prolate spheroid with r — >■ oo, i.e. an infinite rod, we have u zz = 0, and hence 
h = g > 0. For a proof, consider the special case of a dipole density that is fully polarized in 
z-direction with 5=1. The strain be equivalently induced by a set of surface forces 

fi = PijUj. For a cylinder, fi = as n z = 0. Thus, Uij — for S — 1, which implies h — g = 0. 




(5) 
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FIG. 2: The macroscopic force dipole density inside the cellular domain (p^) induces shear strain 
within this domain (v%j) in two different ways, which are characterized by coefficients g and h, see 
eqn. (5). Both coefficients depend on the normalized stiffness E m /E c of the matrix surrounding the 
domain. This is shown for the case of a perfectly spherical cellular domain with r = 1, as well as for the 
case of prolate spheroid with aspect ratio r = 4 and for the case of an infinite cylinder (r = oo). The 
"domain shear compliance" g characterizes the amount of shear strain, which is due to an anisotropy 
of the force dipole density itself. The "shape field factor" h characterizes the shear strain that is 
induced as a consequence of an asymmetric domain shape and anisotropic restoring forces from the 
surrounding matrix by the isotropic "hydrostatic pressure" part of the force dipole density. For the 
figure, the Poisson ratios are v c = v m = 0.3; the limiting cases discussed in the text are marked in red. 

It can be shown that h > holds true for any prolate spheroid with < r < oo. 

Consistent with the limiting cases discussed above, h is a non-monotonic function of E m /E c 
that vanishes if the matrix is either very soft or very stiff. The domain shear compliance g 
is a monotonically decreasing function of matrix rigidity E m since the restoring forces, which 
oppose shear strain w^-, increase with E m . 

The local strain field in the vicinity of a typical dipole at position induced by all the 
other dipoles, u°j C ' k = Yli^k u ij > * s ^ ne superposition of the macroscopic strain field Tiij and a 
local correction term, which is similar to a Lorentz cavity field [9] and insensitive to boundary 
conditions [27]. In a continuum approximation that averages over the random positions x/ of 
the other dipoles [28] , the local strain is induced by a homogeneous dipole density p^ that fills 
the entire cellular domain, except for a spherical void around x& because of hard-core repulsion. 
Thus, the local correction approximately amounts to subtracting the strain contribution of a 
dipole density p^ filling the void and we obtain for the shear strain 

((u l ° c,k,s (x.))) k = m-j - a s plj/E c for |x - x fc | < a, (6) 

By construction, the averaged local strain field looks the same for all force dipoles, except 
possible those which are very close to the boundary of the domain Q. 



7 



r \ Monotonic order-response 

v d / for strong interaction, high density, low temperature 




0.02 0.1 0.2 1 2 10 20 oo 



Non-monotonic order-response 
v / for weak interaction, low density, high temperature 

0-2 I" S 




Oi . I i . > . 

0.02 0.1 0.2 1 2 10 20 oo 



FIG. 3: Two regimes of rigidity sensing. The nematic order parameter S of the force dipole ensemble 
depends strongly on the normalized stiffness E m /E c of the matrix. Dependent on the effective strength 
J/T* of their elastic interaction, we distinguish two regimes for this dependence: (a) We observe a 
monotonic, step-like dependence on E m /E c for strong interactions with J>T*, i.e. for strong dipole 
strength Pq, high dipole density p and low effective temperature T*. A similar dependence is found 
for S as a function of J/T* if E m is fixed (here E m = 2E C ), see inset, (b) We observe a non-monotonic 
dependence with a maximum of nematic order at an optimal matrix stiffness for weak interactions 
with J <C T*, i.e. for weak dipole strength Pn, low dipole density p and high effective temperature 
T*. The Poisson ratios are always v c = u m = 0.3, except for panel (b), where we show in gray also 
the case u c = 0.45. 

V. MEAN FIELD THEORY FOR NEMATIC ORDER 

The Hamiltonian H = ^2 k iUki of pairwise elastic interactions of the force dipoles can be 
rewritten as a coupling between their dipole moments and the respective local strain fields (plus 
a self-energy term that is independent of the dipole orientations) 

H = \Y,I ^^ k P^/\B\ + U kk . ( 7) 

Z k ^|x-x fc |<a 

We derive a mean-field Hamiltonian Hq of Lorentz- Weiss type from this Hamiltonian by replac- 
ing the true local strain fields u° c ' k by their mean field averages ((v 1 ° c ' k )) k . Using eq. (5) (and 
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a standard variational technique, which results in a prefactor 2 [18]), we find 

Ho = J2 [ d^((v}^ k )) k P^/\B\ 

k ^|x-x fc |<a ^ 

= const. — ^ J (h — gS) Sk 

k 

with g = g — u s - Here J = 2pP$ /(3E C ) sets the energy scale of elastic interactions. This 
energy scale should be compared to an effective temperature T* that quantifies the noise in 
the system [19]. Assuming for simplicity Boltzmann statistics for the dipole orientations, we 
self-consistently solve for the nematic order parameter S [20]. In the limit of low interaction 
strength \g\J <T*, we find S « hJ / (5T* + gj). For prolate spheroids, the shape field with h 
is always positive and favors nematic order. The coefficient g = g — a s can be either positive 
or negative as the domain shear compliance g > and the cavity effect of isotropic hard-core 
repulsion (a s ) give competing contributions to g. From eqn. (6), this can be understood in 
terms of increased or decreased structural interference of the strain induced by the "central" 
dipole at with the macroscopic shear strain and the local cavity correction —a s plj/E c , 
respectively. 

The elastic interaction energy of eq. (8) is formally equivalent to the generic theory of 
nematic liquid crystals of Maier and Saupe [21] with an additional, external alignment field 
[20]. In our case, this field arises from the contractility of the system (which has no analogue 
in the electrostatic case) and does not require an external stress. The field strength hJ is 
proportional to the shape field factor h and thus depends on the shape of the cellular domain 
as well as the rigidity of the matrix. Recall that for h = (i.e. for spherical domains or very stiff 
or soft matrices), Maier-Saupe theory predicts a first order phase transition from an isotropic 
to a nematic phase when gJ/T* —4.542 [21]. A necessary condition for this transition is 
g = g — a s < 0, i.e. the cavity effect of isotropic hard-core repulsion has to outweigh the 
domain shear compliance g. In the special case of a spherical domain with r = 1, the shape 
field vanishes, h = 0, and we have g < a s for stiff matrices with E m > E c , see fig. 2. Accordingly, 
we find strong nematic order for E m > E c , see fig. 3(a). The full phase diagram exhibits a 
line of first order phase transition that ends in a critical point, see [20] and fig. 4. Generally, 
positive values of h increase nematic order and smooth out the sharp phase transition that 
exists when h = 0, see fig. 3(a) for the case r = 4. This monotonic, step-like dependence of the 
nematic order parameter S on matrix rigidity E m is found whenever the elastic interactions are 
sufficiently strong with a s J/T* > 4.542. For weak elastic interactions (or large noise strength), 
the nematic order parameter S emulates the behavior of the shape field factor h, and hence 
depends non-monotonically on matrix rigidity E m for prolate spheroids, see figs. 2 and 3(b). If 
the cellular domain Q is a perfect sphere with r = 1 and elastic interactions are weak, there is 
no nematic order. The results shown are robust with respect to changes of the Poisson ratio u m 
of the the matrix, and, in the regime of a monotonic nematic response, also to changes in the 
Poisson ratio v c of the cellular domain (not shown). Only the non- monotonic nematic response 
is attenuated in nearly incompressible cellular domains and vanishes for v c = 1/2, see fig. 3(b). 
In conclusion, the nematic order of force generators provides a read-out of matrix rigidity and 
exhibits qualitatively different regimes of dependency for strong and weak elastic interactions. 
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FIG. 4: Phase diagram for the nematic order parameter S of an ensemble of force dipoles as a function 
of the coefficients of the Hamiltonian Hq, see eqn. (8). In our cell model, the coefficients g and h depend 
in turn on domain shape and matrix stiffness in a non-trivial way, see fig. 2. As illustration of this 
dependence, we traced out these coefficients for an example parameter set for a prolate domain (r = 4, 
v c = v m = 0.3), which correspond to the black E m -S curve in fig. 3(b). 

It should be noted that our mean field theory is based on a Weiss approximation and does 
not correct for a reaction field [22]; we thus underestimate g and overestimate the transition 
temperature. In a real cell, anisotropic hard-core repulsion of actin bundles can favor nematic 
order [23] and thus effectively decrease g. Both corrections stem from local effects and are 
insensitive to boundary conditions. 

A non-monotonic dependence of cell force polarization on matrix rigidity was already pre- 
dicted in [4] in the framework of a linear response theory, which assumed pf» = — Xs^lj- This 
linear response theory corresponds to the regime of weak elastic interaction J <C T* in our 
theory. Our theory allows a self-consistent derivation of the phenomenological parameter Xs, 
which represents an orientational polarization susceptibility: In the limit J < T*, we find 
Xs = {E C /5)J/T* . Other work studied nematic order by elastic interactions within a two- 
dimensional ensemble of force dipoles using Monte-Carlo simulations and found nematic order 
for high dipole densities [12]. The use of periodic boundary conditions in that study is equivalent 
to a strain-free boundary and should be compared to our limit E m 3> E c . 



VI. CONCLUSION 

We showed that elastic interactions can drive nematic ordering of cytoskeletal force gener- 
ators such as short, nascent stress fibers as a function of matrix rigidity. A non-monotonic 
dependence of actin cytoskeleton polarity on matrix rigidity has indeed been found in experi- 
ments with embryonic stem cells [4], which corresponds to our limit of a low density of force 
generators and large noise strength. Interestingly, in high aspect ratio cells, nematic order sat- 
urates at high matrix stiffness [4]; this may correspond to our prediction of saturation at high 
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values of J/T* since these cells display more developed stress fibers and hence higher contractile 
forces. We speculate that different cells might employ different mechanisms of rigidity sensing 
by controlling e.g. the strength of active cell forces. 
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